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SUMMARY 

INTERP3 is a computer routine designed to linearly interpolate a variable 
which is a function of three independent parameters. It is unique in that the 
variables witl.in the parameter arrays do not have to be distinct, or equally 
spaced, and the array variables can be in increasing or decreasing order. 


INTRODUCTION 

The type of tri variate functions which can be evaluated by existing interpola- 
tion routines is severely limited by restrictions imposed by the routines. In 
the Langley computer complex, the.e routines include ITRI, which requires the 
parameter variables to be distinct and in algebraically increasing order, and 
FUN6EN which assumes that the function is known for equal intervals of the 
parameter variables. 

The desire to use the computer for linear interpolation in a data table of 
propeller slipstream characteristics generated the need for a computer routine 
to interpolate a trivariate funccion defined by nondistinct, unequally spaced 
parameter variables. INTERP3 is designed to fulfill that need. Furthermore, 
this routine allows the parameter variables to be in increasing or decreasing 
order. 


Description of Program 

This proaram performs linear interpolation for a set of functions defined over 
the same set of independent parameter points. By the use of multiple entry 
statements, INTERP3 evaluates the ^unction at the point (XIND, YIND, ZIND), 
producing a value WD, where: 



XIND = The input x parameter where interpolation is desired 

YIND = The input y parameter where interpolation is desired 

ZINO = The input z parameter where interpolation is desired 

WO = The output interpolated value of the function at the 

point (KIND. YIND, ZIND). 


This program does not allow extrapolation. If any input parameter exceeds the 
data table range, a warning is printed and WD will not be calculated, but 
given the value of zero. 

The names and descriptions of the subroutines are as follows: 

INTERP3 - Outnn\y main program to initiate calculations. Control 
returns to this routine after evaluation of each input 
point. 

XSEAR - Main routine. Loads data table and performs inter- 
polation in X parameter array. 

YZSEAR - Performs interpolation in Y and Z parameter arrays. 

WSEAR - Performs interpolation in W array. 


A FORTRAN listing of this program is presented in Appendix A, and a flow chart 
is presented in figure 1. 


Description of Data File 

The data file contains all tabular values in the X, Y, Z, and W arrays. The 
purpose of FORMAT 601 is to read these parameter values in the sequence x, y, 
z, w for each record. To ensure that all parameter values are read into the 
program, a record must be added to the end of the data table file with dummy 
values that do not equal the respective x, y, z, w values of the last true 
record. This edited file must be placed at TAPES, prior to execution of 
the program. 

A sample case including data file and input/output listings is presented in 
Appendix B. 


CONCLUSION 

INTERP3 is a program designed to linearly interpolate a variable which is a 
function of three independent parameters. It is unique in that values within 
the parameter arrays do not have to be distinct, or equally spaced and array 
values can be either algebraically increasing or decreasing. These features 


2 



of INTERP3 give it a capacity of interpolating many different kinds of tri 
variate functions. 
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APPENDIX A 

LISTING OF PROGRAM INTERP3 
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APPENDIX B - SAMPLE CASE 


(a) Data File Listing 

The following data file is listed as it would appear at TAPES. 
FORMAT 601 reads the values x,y,z,w, respectively. Note the 
arrangement of the values in each paramter array, and the last 
file record, which contains dumn\y values. 


X 

y 

z 

W 

• 5000 

0.0000 

.2070 

• 3280 

.5000 

o.cooo 

.3420 

.4361 

.5000 

o.cooo 

.5037 

• 5142 

.5000 

0.0000 

.6846 

• 5524 

.5000 

• 2000 

• 2000 

.2905 

.5000 

.2000 

.4000 

.4580 

• 5000 

.2000 

• 6000 

.5660 

.5000 

.2000 

.8000 

• 6140 

.5000 

.2000 

1.0000 

• 6180 

.5500 

3.C000 

.4 863 

.1345 

.5500 

3.0000 

.6777 

.1940 

.5500 

3.0000 

1.0675 

• 3056 

.5500 

3 .0000 

1.4706 

• 4082 

.5500 

3.0000 

1.8861 

• 5030 

.5500 

3.0000 

2.3118 

.5899 

.5500 

3.2000 

. 5994 

.1564 

.5500 

3.2000 

1.0255 

.2756 

.5500 

3.2000 

1. 4662 

.3848 

.5500 

3.2000 

1.9204 

• 4 862 

• 5500 

3 .2000 

2. 3846 

.5766 

• 5500 

3. <.000 

.3006 

.0555 

.5500 

3. <.000 

. 7581 

.1877 

.5500 

3.^000 

1 .2318 

• 3087 

.5500 

3 .AOOO 

1.7204 

.4201 

.5500 

3.<»000 

2.2217 

.5217 

.6500 

3.2000 

• 6146 

.1607 

.6500 

3.2000 

.8316 

.2227 

.6500 

i.2000 

1.0530 

.2825 

.6500 

3.2000 

1 .2782 

.3395 

.6500 

3.2000 

1.5C70 

.3941 

.6500 

3.2000 

1. 7397 

.4469 

.6500 

3 .2000 

1.9751 

.4970 

.6500 

3.200C 

2.2132 

.5441 

.6500 

3.8000 

.4215 

.0719 

.6500 

3.8000 

.9895 

.2160 

.6500 

3 .6000 

1.5758 

.3471 

.6500 

3 .8000 

2.1778 

.4655 

.6500 

3.8000 

2. 7915 

.5724 

0.0000 

0.0000 

O.OOOO 

0.0000 



APPENDIX B - Concluded 


(b) Input/Output Listing 

Compare the following input/output values to those presented in 
the sample data file. Error messages will be printed whenever 
extrapolation is attempted. 


XIN0» 

• 5500 

YINO* 

3.0000 

zin: ' 

1.4700 

WO* .4080 

XINO- 

.5000 

YINO* 

. 2000 

UNO* 

.3000 

WO* .3743 

XIND« 

.5500 

YINO* 

3.3000 

ZIND* 

2.0000 

WO* .5017 

INPUT VALUE 

OF XIND 

1 OUT OF 

X TABLE 

RANGE 


XIN0« 

.7000 

YINO* 

2.8000 

ZINO* 

• 2500 

WO-0.0000 

XINO* 

.6000 

YINO* 

3.2000 

ZINO* 

1 .5000 

WO* . 39 24 

XIND* 

.6500 

YINO* 

3.8C0C 

ZIND* 

1.0000 

WO* .2183 

XINO* 

.6500 

YINO* 

3.6000 

ZINO* 

2.0888 

WD* .5195 
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